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Abstract 

We present a general procedure for incorporating higher-order informa- 
tion into the scale-setting prescription of Brodsky, Lepage and Mackenzie. 
In particular, we show how to apply this prescription when the leading 
coefficient or coefficients in a series in the strong coupling as are anoma- 
lously small and the original prescription can give an unphysical scale. 
We give a general method for computing an optimum scale numerically, 
within dimensional regularization, and in cases when the coefficients of a 
series are known. We apply it to the heavy quark mass and energy renor- 
malization in lattice NRQCD, and to a variety of known series. Among 
the latter, we find significant corrections to the scales for the ratio of 
e^e~ to hadrons over muons, the ratio of the quark pole to MS mass, the 
semi-leptonic B-meson decay width, and the top decay width. Scales for 
the latter two decay widths, expressed in terms of MS masses, increase by 
factors of five and thirteen, respectively, substantially reducing the size of 
radiative corrections. 
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1 Introduction 



QCD processes computed to a finite order in perturbation theory depend on 
both the choice of renormaUzation scheme and the scale for the running cou- 
phng constant asiq). In particular, changes in the scale induce variations at 
the first neglected order. While these variations diminish as higher orders are 
included, for low-order calculations they can be significant, particularly for pro- 
cesses sensitive to relatively low scales. Finding an optimum, physically mo- 
tivated method for choosing this scale in such cases is important not only to 
produce accurate results, but also to reasonably estimate convergence based on 
the size of the series terms. Such a method allows a meaningful prediction or 
comparison with data even at leading order. 

A variety of procedures have been proposed to select this scale [|l|-(l7|. In 
this paper, we investigate the prescription of Brodsky, Lepage and Mackenzie 
(BLM) jj]. In this method, one chooses the scale q* for as{q*) which approx- 
imates the use of the fully dressed gluon propagator within that process. The 
choice is equivalent to determining the dominant momentum flowing through 
the propagator within a diagram It has been applied successfully in 

a large variety of perturbative calculations. Among these, it was essential in 
demonstrating the viability of lattice perturbation theory [ p^ , and in extract- 
ing a precise value of from lattice simulations of the T and V' systems |l^ . 

In this paper, we generalize the prescription to remedy an anomaly observed 
in a variety of applications, particularly apparent when determining the scale 
over a range of parameters in the action. The NRQCD mass and energy renor- 
malizations presented in Sect. |9.2| are typical examples. In most of these cases, 
for some value of the bare quark mass, the BLM scale diverges. We show that 
this breakdown is not a flaw in the general prescription, but rather the result 
of employing only a single vacuum-polarization insertion to estimate the typical 
momentum. While we focus on setting the scale for one-loop diagrams, we use 
information from two-loop and higher insertions within these diagrams to pro- 
vide a simple generalization which accurately estimates the scale over the full 
range of parameters. It is straightforward to implement for both analytic and 
numerical computations, requiring only a modest extension beyond the leading 
order determination. For both computations, one obtains the additional infor- 
mation required from one higher moment in \og{q^) within the same diagram 
as was used in the lowest order application. For processes where the series 
coefficients are known, it requires only identifying the coefficient from vacuum 
polarization at the next order. 

We note that other authors have developed a variety of extensions to Ref . Q , 
which explore conformal symmetry and the relation between various perturba- 
tive schemes ^, 0, or which estimate nonperturbative contributions and 
resum classes of diagrams to all orders [^4|, 16 . Our goal is more modest: to 



provide a simple but robust scale determination for a process calculated to finite 
order. Specifically, we choose a single optimized scale for the leading, one-loop 
diagram, to be used for all orders. We show, however, that our prescription 
should effectively absorb into the leading term or terms the bulk of contribu- 
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tions from all higher order diagrams which dress the leading gluon. 



2 General prescription 

Following Refs. § we choose the V scheme based on the static-quark po- 
tential because of the direct connection between the scale of its coupling ay and 
the momentum flowing through its associated gluon. For a one-loop diagram 
with an integrand f{q) which contributes predominantly at large q, a natural 
choice for the scale q* of av is a mean value which reproduces the result of a 
fully dressed gluon within the diagram j|, |l^ , 

aAq*) I d\fiq) = [ d\ a,.{q)f{q) , (1) 



as illustrated in Fig. |^. However, av{q) possesses a pole at Ay, an artifact of 
an all-orders summation of perturbative logarithms. We avoid this singularity 
by truncating the series for av{q) at a finite order, as is appropriate for an 
asymptotic series. 

Expanding av{q) in terms of ayiq*) 

= 1^ r^;i'^W2/.2^ - - aliq*)f3ologiqVq*') + ■ ■■ (2) 

and solving to first nontrivial order gives ||l3|] , 

, /A/(g)iog(g^) (/(g)iog(g^)) //, , 2,. 
log(<z ) = ;,^;(,) - (/(,)) ^ )> ' 

a statement of this prescription suited for numerical calculations. Here 

and indicates an average weighted by /(<?). 

By the definition of a^, Eq. (||) guarantees that av{q*) absorbs the effect 
of second-order vacuum polarization insertions in the gluon's propagator. An 
alternate method to determine q* is then to require that [3q, or equivalently Uf, 
disappears to that order Q . This version is useful when the /3o orn/ dependence 
of coefficients in a perturbative expansion are known explicitly. We discuss this 
in more detail in Sect. ^. 

Eq. (H) produces an optimum scale by means of an average of log((7^) weighted 
by f{q). As such, it provides a measure of the typical momentum carried by this 
gluon in the dominant integration region, in accord with intuition. However, 
in certain cases (/) vanishes, rendering q* from Eq. (||) meaningless. This is 
a consequence of using an expression first-order in av{q*) for a process which 
is properly second-order, rather than a flaw in the general prescription. When 
(/) vanishes, the diagram on the left in Fig. |l| does not contribute. The leading 
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contribution from this gluon is second order, and the requirement that q* be 
chosen to best approximate the all-order result leads to the equation illustrated 
in Fig. 1^. For an integrand dominated by large momentum, the left side of 
Eq. (pi) is replaced by 



-aUq*)Po J d\f{q)\0g{q'lq*') , (5) 

as is known from the running of av{q). Expanding av{q) as in Eq. (|^) yields 

(/log^(g^)) 

This, rather than Eq. (|^), is the appropriate statement of the prescription for 
this case. 

As a simple illustration, consider a model in one dimension in which the 
Feynman diagram produces an integrand 

f{q) = 5{q - qa) - 6iq - q^) , (7) 

with positive qa and qt- A reasonable expectation would be that q* should 
be some average of the contributing scales qa and qt, particularly if they are 
nearby. Because (/) vanishes identically, Eq. (||) produces a divergent q* , 
whereas Eq. (|6|), which begins with the next order contribution, gives for q* 
the more reasonable geometric mean 

q* = Vlaqb ■ (8) 

This is the same scale obtained by Eq. (^) applied to the positive integrand 

f{q)^Siq-qa) + Siq-qb) , (9) 

as might be expected. 

In other cases, while not strictly vanishing, (/) may be anomalously small, 
and the dominant contribution from this gluon is still as in Fig. ^j. It is useful to 
generalize the lowest order prescription of Eq. (|^) to incorporate both these cases 
naturally, and also to anticipate the situation where (/log(q^)) is anomalously 
small. 

The discrepancy between the left and right sides of Eq. (|l|) relative to the 
lowest order term av{q*) J d\f{q) is 

- a.(<7*)/3o(log(gV'7*'))) + aliq*)P'oi^og'{qyq*')} (10) 
Applying Eq. (^ leaves a leading difference of 

aliqimogiq') i\og{q')}]'} ^ al{q*)PW (H) 

with cr the standard deviation of log(g^) with respect to the weight f{q)/ (/). 
When f{q) does not change sign, f{q)/ (/) and are both positive. When f{q) 
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Figure 1: The BLM prescription for fixing the optimum scale q* to leading 
order in as{q*). 




Figure 2: The BLM prescription applied to a process for which a gluon con- 
tributes first at order a^(q*). The insertion on the left side represents vacuum 
polarization from both quarks and gluons. 



changes sign and (/) is anomalously small due to cancellations, this error can 
become arbitrarily large, indicating that treating this as a first order process is 
invalid and it is useful to incorporate information from the next order. 

Matching the gluon's contribution to next order to the fully dressed gluon, 
as in Fig. ||, adds the term in Eq. (||) to the left side of Eq. (^. Expanding both 
sides in terms of av{q*) as before leads to a leading relative difference of 

al{q*)PU{^og\q'/q*')} = (12) 
-al{q*)f3'o [\og\q*') - 2(l0g(<Z^))) \og{q*') + (log^iq')}] . 

When f{q)/ (/) is positive for all q, this discrepancy is also strictly positive, 
and the best that can be done is to choose q* to minimize it. The result of 
minimization is again just Eq. (^). This will also clearly be the case when /(g) 
changes sign in some small region without significant cancellations. The leading 
error is then the same as in Eq. (pi]). 

However, when f{q) possesses significant sign changes, the error in Eq. ( p^ 
can become negative for certain values of \og{q*^), and minimization is not 
appropriate. In this case, it is possible to eliminate the difference altogether by 



5 




Figure 3: The BLM prescription applied to second order in as{q*). The 
insertion on the left side represents vacuum polarization from both quarks and 
gluons. 



choosing one of the two solutions 

(/log(g2)) ± [(/log(g2))2 _ (/) (flog'iq^)) 



(/) 

^ (log(g^)))±M^ (13) 

When the logarithmic moments are available over a range of parameters, re- 
quiring that q* be continuous and physically sensible makes the proper choice 
apparent, as will be observed below. In particular, when (/) is nearly zero, this 
requires choosing the sign opposite to that of (log(g^)). In every case we have 
considered, the choice has been obvious. However, if the need arises, one may 
resolve the sign unambiguously by using information from higher moments, as 
discussed in Sect. 0. 

Were f{q) a probability distribution, a would be its standard deviation. A 
negative value for cr^ indicates that f{q)/ (f) has substantial changes of sign 
and is behaving significantly unlike a probability distribution. In this case, (/) 
will be anomalously small, the order a'^{q*) contribution becomes important, 
and Eq. (|l^) provides the appropriate choice of scale. As a result, Eq. (|3|) 
determines q* when cr^ is negative, Eq. ^ when positive. This prescription is 
the main result of this paper. 



Although Eq. (13) uses information from part of the next order, it is the 
appropriate scale to use when is negative, even if only computing to first 
order in av{q*)- In that case, in addition to setting the scale for the leading 
term, it allows for a reasonable estimate of the magnitude of the neglected next- 
order terms, based on a^(g*). When (/) is very small, these neglected terms 
should give a sizable correction to the first-order term. And when one computes 
to order a'^{q*) or higher using this scale, higher-order terms which dress the 
leading gluon should be small, having been largely absorbed into the first two, 
as in Fig. p|. 
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Figure 4: The BLM scale q* for the model of Eq. ( p^ as a function of c, with 
qa = 2.0 and qb — 1.8. The first order solution of Eq. ^ determines q* for 
c < — 1, the second order solution of Eq. ( p^ for c > — 1. The dark dotted lines 
show the first-order solution in regions in which it does not apply; light dotted 
lines display inapplicable second-order solutions. 

As an illustration, we consider a slightly more general version of the model 
OfEq. (0), 

f{q) = {l + c)S{q-qa)-S{q-qt) , (14) 

where (J) vanishes exactly for c — 0, and has partial cancellation for c > — 1. 
Figure ^presents the scale q* determined by Eq. (||) when cr^ > and by Eq. (O) 
when (T^ < 0; that is, when c > —1. For c < —1, with no cancellations, Eq. (^) 
produces reasonable values for q*. It falls between qa and qb, approaching qa 
for |c| large, and q^ for c = — 1. As c approaches zero and (/) vanishes, q* from 
this prescription diverges. However, for c > —1, cr^ is negative, and Eq. (^3|) 
provides the optimum scale. It evidently behaves according to expectations. 
Even for the case when c is large and positive and Eq. (^ produces a fairly 
sensible result, it overestimates q* and Eq. ( p^ is preferable. 

To summarize, we restate the prescription in a more compact form. We have 
chosen inclusion of the running coupling within the integrand for a first-order 
diagram, / d\ av{q) fiq), as a natural means to account for the running of the 
coupling with the gluon's momentum. It has the advantage that it incorporates 
higher-order diagrams which dress the gluon, has no arbitrary scale dependence, 
and appropriately accounts for the strength of the coupling of a gluon with 
momentum q. 

It has the disadvantage that av{q) has an unphysical pole at q = Ay, making 
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the integrand ill-defined. We avoid this by expanding av{q) at the scale q* as 
in Eq. (|^), and working to finite order in av{q*), giving 

J d\f{q) K(g*)-/3o4(9*)log(gV9*')+/3o«v('?*)log'(<zV9*') + ---] • 

(15) 

We choose q* to reproduce the full integral as well as possible. In the absence 
of significant cancellations in / d\ f{q), we may select the scale by Eq. so 
that the first non-leading term in Eq. (^5|) vanishes. The discrepancy is then 
the term of order a^(g*), which this choice for q* minimizes. Furthermore, as 
q* will be near the typical q, f{q) will be roughly even about q* , and higher- 
order contributions should be either near zero or their minimum depending on 
whether they are even or odd in av{q*)- 

However, when f{q) is essentially odd about some q and so suffers from 
significant cancellations, this is not an appropriate prescription. The leading 
term in Eq. (|l5|) will be anomalously small compared to the second term; in 
extreme cases, it might even vanish, and it would no longer make sense to absorb 
the second term into the leading term. Furthermore, the scale from Eq. (||) 
would no longer accurately represent the typical momentum, and neglected 
higher order terms in Eq. (|l^) would be anomalously large. It is, however, 
possible and sensible to require the third term to vanish by Eq. (^; that is, to 
absorb it into the second. This again provides a typical scale about which, in this 
case, f{q) is essentially odd, minimizes the fourth-order term, and suppresses 
higher-order terms. 



3 Schemes other than V 

For prescriptions other than a^, vacuum polarization insertions will in general 
contribute subleading constants in addition to terms as in Eq. (|). Though 
nonleading, these constants can make significant contributions at physically 
interesting values of g^, and so the optimum scale ought to be chosen to account 
for them as well. One method for doing so is to focus on the fermion loop Q. 
Both the log(g^) and the subleading constant will appear multiplied by iif. 
Replacing n / with /3o using Eq. (^) modifies the fermion loop contribution by a 
constant, to 

-a,(g*)/3o(log(<zV'Z*')+«) ■ (16) 
When applying the first-order prescription, amending Eq. (|^) to 

\0g{q*^) = {\0g{q^) + a} (17) 

absorbs both the leading log and subleading constant into as (<?*)■ For MS, 
a = -5/3 resulting in the shift in scale g 

qt^ = exp(~5/6) ql = 0.43 ql . (18) 

As with Uv , this also absorbs the log associated with giuon vacuum polariza- 
tion, since /3o determines its contribution relative to the fermion loop. However, 
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Figure 5: One loop self-energy diagram in the g4>'^ model. 



the gluonic subleading constant need not contribute in this ratio, and so will not 
also be completely absorbed. One might choose instead to completely absorb 
the gluonic constant by solving Eq. (|^) for the adjoint Casimir constant Ca = N 
associated with the gluon loop in terms of /3o, before absorbing the /3o term into 
as{q*). This would be particularly appropriate when = 0, for example. For 
MS, the result is a factor of exp(— 31/66) = 0.63, not greatly different from 
Eq. (|l8|). This indicates that to one loop, absorbing the fermion loop constant 
also largely accounts for the gluonic constant. 

When applying the second-order prescription, a constant subleading contri- 
bution leads to the same shift in Eq. ([L3h as in Eq. @, with 



The second term on the right is invariant under a shift in \og{q^) by a constant, 
and so remains unaffected. 

Because the V scheme associates gluon exchange with a physical process at 
the specific scale g*^, higher order contributions associated with the running 
of av must vanish identically when the gluon's momentum hits q*^ . As a re- 
sult, logarithmic contributions \og{q^ / q*"^) from these diagrams appear without 
subleading constants. 



In order to provide a more realistic example, and to show how this prescription 
can be applied simply when using dimensional regularization, we determine q* 
for the one-loop g(jr diagram of Fig. H While we use this as a simplified model 
for a quark self-energy diagram, we also note that this scale setting method is 
not restricted to QCD. 

By introducing an additional denominator of the form {q^ Y into this diagram 
in n-dimensional Euclidean space. 



log(g*2) ^ (log(g2) + ± \(\og[q^)f - {\og\q% " . 



(19) 



4 Determining in MS 




(20) 
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and expanding the result in S to second order, we produce the necessary loga- 
rithmic integrals: 



ifiq)) S (fiq) log(g2)) + - (fiq) log'{q^)) 



(21) 



(Note that Refs. Q present a different and elegant technique for extracting 
these logarithmic moments based on a simple dispersion sum over a fictional 
gluon mass.) 

We evaluate Eq. (E3) using standard methods and obtain 



9^ nS + e) 
(47r)2 T{5) 



dxdyy'-\l-y) 



(22) 



where x and y are the usual Fcynman parameters, e = (4 — ri)/2, /i is introduced 
to keep g dimensionless, and 



M 



P 



(1 — j/)to^ + x{l — y){l — X + xy)p^ 



(23) 



The integration region y ~ produces a l/i5 singularity. Partial integration 
of y^~^ makes this explicit, and allows us to expand in 5 under the integral. 
Keeping terms in 1 /e to order zero and comparing to Eq. ( ^l[ ) gives 



dx log 



+ a;(l — x)p^ 

7^ 



(24) 



log(g^)) = f^dx I - ^ + 1 log (^-|) + ^ 



(25) 



-log^(^)+-log^(^ 



a;2(l — y^p^ 
NP 



{f{q)\oe{q-)) = 2jyxdyi^- ^ + ^ log 



(26) 



TT 

12 



log' 



3 



1 - 



>"(1) 



log 



x^{l — y)^p^ 
JP 



y 



f2 



log' — 



In the above, we have dropped the overall factors g'^ / (Att)'^ , which are here 
irrelevant, and substituted jj? /A-k for /i^. The latter greatly simplifies these 
expressions and allows us to apply the MS prescription by subtracting the e 
poles alone, as one would for MS. 

To renormalize these terms, we note that the \og{q^) and \o^{q^) factors 
integrated within this one-loop diagram stand in for the large-momentum con- 
tributions of the first two higher-order vacuum polarization subdiagrams which 
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sum to form the running coupling. In particular, as these factors are finite, 
they represent these subdiagrams with their subdivergences already removed. 
For example, the factor log{q'^), which appears in the integrand which pro- 
duced Eq. (p5|), comes from the large-momentum approximation to the MS- 
renormalized one-loop vacuum polarization diagram; that is, to Eq. (|2^) with- 
out the pole. The poles in e which remain in Eqs. ( |2^ ) to ( p6| ) are then the new 
overall divergences associated with one, two and three loops, respectively. In 
the MS prescription these are simply discarded. 

Finally, we note that Eq. (^) in this model is the one-loop vacuum polar- 
ization diagram in addition to being the particle's self-energy. At large p^, it is 
approximately logijP' / n"^) — 2, including the subleading constant. The constant 
a in Eq. ( p7| ) is then —2, and the MS value for q* will differ by a factor exp(— 1) 
from the expression for q* in the V scheme. 

While Eqs. (|2j) to ( p6|) allow us to determine q* for any p and m, it is 
illuminating to consider two limits. When ^ , after renormalization and 
choosing /i^ 



P 



ifiq)) = -2 



(27) 



(log(g^))) 



log(/) - ^ 



(28) 



(log^(g2))) = log2(p2)_^iog(p2) + 



TT 
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to leading order in m? /p^ . These give 



V>"(1) _ Tf_ 

~ 12 



TT 

1 + — 



(29) 



(30) 



or numerically, cr^ = .6077. The lowest order solution is then appropriate and 
yields a value for q* close to p, with 



q*/p ^ exp(-l) exp(-7rV48) = .2995 , 



(31) 



which corresponds to a value of exp(— 7r^/48) = .8141 in the V scheme. We can 
use cr to give a measure of the relative spread in the momenta which contribute 
to this diagram, with 

Ag/q w cr/2 = .3898 , (32) 

independent of the prescription. In general, when p^ 3> but for /i^ arbitrary, 
Q* ~ ^/P^■ This is the same result as Eq. (^, and reasonable for a diagram 
dominated by momenta between these two scales. 

In the large mass limit m 3> p, with /i^ chosen to equal the natural scale 
m^, we find 



1 p2 

(/('?)> = t:^ 



(33) 
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(/('?) log(9^)) 



In rrt ^ 



(34) 



(/(g)iog'(g')) = 



(35) 



-1 



6 



■log {rr?) 



to order /rn? . Clearly in this limit {f{q)) becomes anomalously small, and 
we expect the second order solution to be necessary. We confirm this by noting 
that to this order 
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18 



P 



(36) 



which is negative in this limit. The second order formula applies, and gives 



log(g*Vm^) 



2V'"(1) 



3(7r2 - 4) 



9(7r2 - 4) 



(37) 



+ 



27(7r2 - 4)3 



(-3(^2-4)+V/'(l)) 



or numerically, 



q*/m = 0.6953 + 0.0574 p^/m"^ . 
For the V scheme, the leading —2 in Eq. (|37| ) is absent, and 

/m = 1.8899 + 0.1562 p^/m'^ . 



(38) 



(39) 



Fig. ^ and Fig. |^ display q* as a function of p for the respective cases where 
^ = p and II — m. The limiting values discussed above are evident. For = p, 
the first-order solution determines q* in both the large and small p regions, 
connected by the second-order solution in the interim. Immediately to the right 
of the point where the first-order solution diverges in Fig. ^, indicated by the 
vertical line, the second-order solution with positive root determines q*; to the 
left, the second-order negative root applies. For fi — m (Fig. the first-order 
solution applies for large p, the negative root second-order solution for small 
p. In both cases, use of the appropriate second order solution where applicable 
gives a meaningful and continuous value for q* over the entire region in p. Which 
second-order solution to choose from Eq. ([l3| ) is obvious. 

Finally, we note that computing higher order average logs for this diagram 
requires only expanding Eq. (^2|) to higher orders in S, without the need to 
compute additional diagrams. 
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p/m 

Figure 6: The MS BLM scale q* /p as a function of momentum p/m for the 
diagram of Fig. |^ in the scalar (j)^ model, with ^ = p. The first-order solution 
determines q* in both the large and small p regions, connected by the second- 
order solution in the interim. The vertical line indicates the point where the 
first-order solution diverges. The dark dotted lines show the first-order solution 
in the region in which it does not apply; light dotted lines display inapplicable 
second-order solutions. 
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Figure 7: The MS BLM scale q* /m as a function of momentum p/m for the 
diagram of Fig. ^ in the scalar (fi^ model, with /i = m. The first-order solution 
determines q* in the large p region, while the negative root second-order solution 
gives q* for small p. The dark dotted line shows the first-order solution in the 
region in which it docs not apply; the light dotted line displays the inapplicable 
second-order solution. 
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5 Determining q* from a known series 



To apply this prescription we need the first two logarithmic moments within 
the integrand associated with a gluon's propagator. Under certain conditions, 
we may apply this prescription to set the scale for a process for which the 
expansion is already known by examining its Uf dependence. At each order of 
av{lj), the contribution from vacuum polarization will give the largest power 
of n/ , or equivalently, of /3o 0, • It is therefore possible to read off the 
logarithmic integrals directly from the series coefficients. 

Using Eq. (Q) to replace the largest-n/ terms with f3o in contributions asso- 
ciated with a particular gluon, we obtain a series of the form 

coavip) + (ai - ci/3o)av(M) + (^2 H h C2/3o)q^(^) H . (40) 

The coefficients c„ are then associated with vacuum insertions in the gluon 
propagator. Comparison with the right sides of Eqs. (||) and (||) gives 

CO = (/) (41) 
ci/co « (log(gV/^'))) 

which holds when f{q) contributes predominantly at large q. Given this associ- 
ation, the prescription to second order is 

log(g*V/^') = ci/co ± [(ci/co)^ - C2/C0] (42) 
when the argument of the square root is positive, and 

log(g*V/x2) = ci/co (43) 

otherwise. 

For schemes other than V, the presence of a subleading constant a contri- 
bution to fermion vacuum polarization leads to the identification 

ci/co « (log(qVA^2)+a)) (44) 

C2/C0 « ((log(gVA*') +«)'))• 

Because ci includes a and the square root is insensitive to it, Eq. ( p3| ) and 
Eq. ( p2[ ) automatically incorporate the shift in Eq. (|l^) and so may be used 
unchanged. Also, if one is able to identify in the series the constants Ca associ- 
ated with gluonic vacuum polarization, one could choose to use this instead to 



rewrite the series in terms of /3o. Eq. (43) and Eq. (42) would then automatically 



absorb the subleading gluonic constant, as discussed at the end of Sect. g. 

6 Combining series 

Determining the scale for the series formed by multiplying two series, 

Fa^l+Caa^{ql) + --- (45) 
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and 

Fb^l + Cba,.{q;) + --- (46) 

with known scales is straightforward when considering only first-order scale 
setting: 

Fab = FaFb = l + {ca + Cb)a^iq:b) + ' ' ' (47) 

with 

log(.:g)^ ^°^°^^'^"^:^-^"^^'^"^ (48) 

Ca ~r Cb 

Because of the need to first test the sign of the new ct^j^, the prescription for 
applying second-order scale setting is slightly more involved. In that case, 



, _ c,(log^(g^)))^ + c,(log^(g^))), ^ c,(log(g^)))^ + c,(log(g^))), V 



caal + cbal , CaCb .//i,g(^2)^^^ _ (log(<z^))),)^ . (49) 



2 



Ca + Cb (Co -I- Cb) 

As usual, if (T^j > 0, the first order combination 



iog(«:^) = '"'^ ■ (50) 



Ca(log(g2)))„+^(log(g2))) 

Cb 

appHes. If fj^j, < 0, 



iog(gah) ; ± [-'^ab\ ■ (51) 

Ca+Cb 

When combining two series by division, F^^/b = Fa/ Ft, the first-order coeffi- 
cients subtract rather than add. The above formulas again apply, but with the 
replacement Cb — Cf,. These also apply to series combined by addition and 
subtraction, respectively, because the results at first order are equivalent. 

For schemes other than V, one should amend the average logs to include the 
subleading constants, as discussed in Sect. ^. The relation between the scale in 
V and in other schemes remains the same. 



7 Higher orders 

Extending this prescription beyond second order is relatively straightforward, 
though it requires computation of ||log"^(q^)| and higher moments, or informa- 
tion from third-order and higher terms in a known series. An extension to third 
order, for example, would be necessary should both (/) and (/log((7^)) vanish, 
making the third term in Eq. (|l^) the leading term. Absorbing the subsequent 
term by requiring ^/ log'^ ((7^/(7*^)) to vanish would give 

lori,-)^!^, (52, 

3 (/log (g2)) 
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as one would also obtain from Fig. g with two loops on the left side. 

When (/) and \og{q^)^ are not identically zero but are anomalously small 
relative to higher moments, requiring ^/ log'^(g^/g*^)) to vanish still gives the 
appropriate scale. A symptom that this is the case would be a third-order scale 
near Eq. ( |5^ which shifts significantly the scale obtained at lower order, and 
which is more in line with physical expectations. If available, higher even mo- 
ments near their minima or odd moments near zero at this scale would confirm 
it. 

An order-n equation is necessary to set q* only when all of the first {n — 2) 
moments vanish or are anomalously small. It would be unusual for a generic 
integrand f{q) to be effectively orthogonal to more than a few powers of log(g^), 
and so the need to use a high-order equation should be rare. We have found no 
realistic cases for which either Eqs. (H) or (|l^) were not sufficient. 

In Fig. H we illustrate the appropriate scales for a model, 

f{q) = AS{q - qa) + BS{q - q^) + {D + c)S{q - qd) , (53) 

with 

B = ^A\og{qaJqd)/\og{qt>/qd), D^-{A + B), (54) 

contrived such that both (/) and (/log((?^)) vanish at c = 0. The second- 
and third-order solutions behave as expected where appropriate; the first-order 
solution, while not divergent, is low throughout. The unphysical behavior of 
the first-order solution, as well as the significant discrepancy between first- and 
third-order scales indicate that the first-order result is inadequate. 

Note the one- and two- node structures of the integrand f{q) in the two- and 
three-delta models of Eq. (|l^ ) and Eq. ( ^ ) , similar to generic first and second 
excited-state wavef unctions. This is the result of choosing /(g) to be orthogonal 
to the zeroth, and to both the zeroth and first powers of log((7^), respectively. 
Integrands requiring higher-order equations would necessarily have more nodes 
and additional detailed structure. 

In general, higher-order solutions can confirm that a scale determined at a 
lower order is indeed typical. For example, in Fig. |^ we include the third order 
solution for the Eq. (|l^) model; where applicable, it does not differ significantly 
from the first-order result. This will also be apparent when we examine higher 
moments for certain processes in Sect. ^ 

Thus far, for simplicity, we have restricted our discussion to contributions 

to 

Q!v (g) from one loop vacuum polarization. We show here that the above 
expressions for q* are not limited to these. Specifically, including subleading 
contributions to av{q) expanded within a diagram as in Eq. ( p^ ) gives 

aA.q*){f)+o?,{q*)[PaA,]^al{q*)[Pl^2 + Pl^l] (55) 
+ aliq*) [/JflAs + ^/3o/?iA2 + /32Ai] 

+ al{q*)[Pt^i + y/3o'/3iA3 + 3/30P2A2 + |/3?A2 + +■■■ . 
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Figure 8: The BLM scale q* for the model of Eq. ( ^3|) as a function of c, with 
Qa = 2.0, qh = 1.8 and qd = 1.6. The second-order solution determines q* for 
c < — 1, the third-order solution for c > — 1. The dark dotted line shows the 
first-order solution; light dotted lines display inapplicable second and third-order 
solutions. 
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Figure 9: The BLM scale for the (5- function model of Eq. (|lj) as in Fig. ^, 
but with the inclusion of the third-order solution. It is applicable for c < — 1 
and indicated by the additional dark line. Dotted lines indicate inapplicable 
solutions. 



Here 
and 

1 /2857 5033 325 A 

and wc have defined the moments 

A„ ^ if log- {q*yq')) = (/ [log(g*2) _ log(g2)]"^ . (58) 

We see from Eq. ( |55| ) that the reasoning which led to the prescriptions in 
Eq. (HI) and Eq. ([l^) was not dependent on the leading-log approximation to 
aviq)- the diagrams in each new set introduced by increasing the order of the 
leading-log approximation have as coefficients the same sequence of moments 
A„ as in the leading approximation. For non-anomalous cases, requiring Ai to 
vanish via Eq. (^) absorbs the first leading- log correction, proportional to /?oAi, 
into the first term. It also absorbs the first next-to-leading log contribution, 
proportional to /3i Ai, and so on for each logarithmic order. And in each order, it 
leaves a leading discrepancy proportional to the moment A2 , which is minimized 
by Eq. (^. Insofar as this q* represents the typical momentum carried by this 
gluon, higher-order terms proportional to higher moments should be suppressed. 
As a result, the first term av{q*) (/), with the choice of a single scale for q* , 
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should do reasonably well at approximating the right side of Eq. (|]) regardless 
of the number of loops kept in the (3 function for av{q)- When (/) vanishes or is 
anomalously small, it is inappropriate to absorb the leading correction at each 
logarithmic order into the vanishing or small first term. In this case Eq. ( p^ 
gives the appropriate scale. It causes A2 to vanish, and so absorbs the second 
correction at each logarithmic order into the set of leading corrections. Again, 
because Eq. (^3|) produces a typical scale, higher-order terms should be small. 

This illustrates one of the advantages of using a coupling based on a physical 
process, such as av For other schemes, there can be significant subleading 
constants associated with the diagrams which dress the gluon, not accounted 
for in the running coupling. These will appear, for example, in the coefficients 
associated with the leading log term at each order n; that is, with /3o ~^, as 

A„ ^ (/ [log(g*2) _ (iog(q2) + ^ (59) 

and one may absorb them by adjusting the scale. However, only one such 
constant can be absorbed in this manner; parts of those associated with j3i and 
higher will remain. By its definition, these constants cannot appear for a^, and 
av{q*) well represents the strength of a physical gluon at the scale q* . This 
optimum choice of scale minimizes all coefficients associated with dressing the 
gluon, not just those associated with j3o. 

8 Improving convergence 

Thus far we have been considering the appropriate scale q* for the leading ay. 
For the sake of simplicity, we will be content to optimize the scale for the leading 
term, and will use the same scale for higher-order diagrams. From the previous 
discussion, it is clear that this will be a reasonable scale for diagrams which 
dress the leading-order gluon; these contributions should be small, having been 
largely absorbed into the leading term or terms. There is no reason to believe, 
however, that it will be the best scale for other higher-order diagrams, and it 
should certainly be possible to improve the convergence of the series by choosing 
the scale for such diagrams separately ||, |lO) . 

There are other cases for which it could prove advantageous to allow different 
scales at different orders, for different diagrams within the same order, or even 
within a single diagram. Returning to the simple and unexceptional model of 
Eq. (^), for which the first-order scale setting Eq. gives Eq. (^), we note that 
the moments 

{\0g-{q*^/q^))) = \ [W^cf-^lql) + log"(g*V%')] (60) 

become zero for n odd, and \o^[cial qh) for n even. The latter grow in magnitude 
when the qa is a few times greater than q\,. As these are proportional to the 
coefficients of terms which dress the gluon, terms in the series become « 1 when 
this range qa/qb exceeds a few over ag. The series becomes badly behaved then 
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not only when the scale for is low, but also when the range of important 
momenta is large |l^ . 

Nevertheless, in this simple case the remedy is obvious. The problem results 
from requiring as at a single scale to incorporate two widely different scales. 
Separating these and writing the series for this model in terms of both as{qa) 
and as{qb) incorporates vacuum polarization contributions exactly and causes 
higher moments to vanish. This reinforces the idea that for a series in which 
gluons from different diagrams occur in loops sensitive to significantly different 
momenta, allowing the associated with each to have its own scale could 
improve the series' convergence. Furthermore, for a series in which a gluon in a 
single diagram is sensitive to a wide range of momenta one might even consider 
improving its behavior by splitting up the integrand, with at a different scale 
assigned to each region, as in the example above. 

Even for diagrams which dress a specific gluon, it is possible to minimize 
higher moments by allowing the scale associated with these diagrams to differ 
at different orders. In particular, one may select the scale at every other order 
such that the following moment vanishes, and the next is minimized. Variations 
in the scale would account for the different regions the moments probe in the 
integrand. 

We will not pursue this further here, but mention that care must be taken 
to preserve gauge invariance if separate scales are assigned to different parts of 
a series. 

9 Applications 

9.1 Known series 

In Table |^ we present a collection of results for perturbative quantities for which 
at least the second logarithmic moment is available, allowing us to apply scale 
setting beyond lowest order. Ref. ||l^ presents a useful compilation and discus- 
sion of many of these. These include the log of the 1x1 Wilson loop in lattice 
QCD (— log W^ii) the ratio of e+e~ goes to hadrons over muons (i?e+e- ) 

[26 -pO{ , the ratio of the quark pole mass to its MS mass (M/M) |^-[^ and 
[14, |16| , the ratio of t goes to i^r + hadrons over r goes to Vre'T^i (Rt) [0~ 
[12 and j2^, |3^, |l^, the scmileptonic i3-meson decay width {T{B — > XueV)) 
[13 -|5^ expressed both in terms of the pole and MS 5-quark masses, the top 
quark decay width (r(< bW)) and [|4[ |6[ [T|] in terms of its 

pole and MS masses, the Bjorken sum rules for polarized electroproduction 
- gf{x,Q^)]) [|5^-[Q, and deeply inelastic neutrino-nucleon 
scattering (/ dx\Fj'^(x, Q^) — Q^)]) [^-[Q, and the static quark po- 

tential iV{Q^)) ||-[|| and 

We note that the non-singlet part of the Ellis- Jaffe sum rule [^-[^ and 
[ |62| gives the same scale as the former Bjorken sum rule, the Gross-Llewellyn 
Smith sum rule [^, ^ the same scale as the latter. 

All but the first scales are from known MS series, and for these at least 
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the fermion vacuum polarization graphs must be given to two loops. In several 
cases higher logarithmic moments are known, which will allow us to test the 
consistency of the procedure. We find four cases where the second-order formula 
gives the preferred scale. 

When the first-order solution is appropriate, the second moment gives a 
rough measure of the range in momenta which flow through the gluon, with 



Aq 1 /A2 
Q "^2 V(/> 



(61) 



Here, Aq is the standard deviation in q and A2 is defined in Eq. (pq). A large 
range results in large coefficients at higher orders, as discussed above. When 



the second-order scale is appropriate, Eq. (61) clearly is not. However, if higher 



moments are available, we may estimate this range using 

(62) 



Aq 

q 



4(n+l)A, 



with n odd. In Table |l|, we use n — I. This expression gives the standard 
deviation in a distribution modelled by a Gaussian times [log(g*^) — log(g^)] 
to render it odd. We found it to give reasonably consistent results for various n 
when applied to several examples discussed below, though other measures are 
certainly possible. 

For R^+f,-, M/M, and both T{B XucU) and r{t — > bW) expressed in 
terms of MS masses, we find that the second-order scale is appropriate, leading 
to significant corrections to the anomalously low first-order scales, especially in 
the latter three. While the new scale for M/M is significantly increased, we 
note that Aq/q is still relatively large, indicating sensitivity to low-momentum 
scales even when AI is large, and threatening a poorly behaved series. This 
apparently infects the b and t decay rates when expressed in terms of pole 
masses, as shown by their low scales. By contrast, MS masses behave more as 
bare masses, being sensitive to short distances; expressing the two decays in 
terms of these significantly improves their behavior js^, [tJ, |2j, |lj, This 
is clear from both their scales and widths. Both these series should be well- 
behaved, and well-represented by as at a single, physically reasonable scale. 
But it is necessary to use second-order scale setting to see this; the first-order 
q* for each indicates a scale which is misleadingly low. 

Refs. 11, |l|, provide very useful values for fermion vacuum polarization 
contributions, and therefore logarithmic moments, computed to eighth order for 
the pole to MS ratio, r, B and t decays. These allow us to compute their Ag's 
using Eq. ( |6^ , but more importantly, to confirm the general picture as discussed 
m Sect. 0. In Fig. ^ we use this information to display the first eight moments 
|( [log(g^l) - log(g2)]"_^^/" as functions of log(g^/M2) for r{B X^eV,) ex- 
pressed in terms of the MS mass Mb- Here q^ absorbs the fermion loop constant 
associated with the MS prescription, as in Eq. (p^), and Mb is the fe-quark pole 
mass. We observe that choosing log(g^/M^) to set the second moment to zero 
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using Eq. ( p^ ) not only removes it and minimizes the third moment, it also sets 
all of the higher moments near their minima or zeros. It is clear that this is 
the natural scale for this process, and that terms beyond second order which 
dress the leading gluon should be small. The first moment is clearly anomalous, 
and setting it to zero using Eq. would evidently lead to large higher-order 
corrections. In general, we expect f{q) to be either roughly even or odd about 
its typical scale q*, and the sign of the second moment, cr^, should distinguish 
the two. For Aq sufficiently small, using Eq. (^) or Eq. ( p^ depending on the 
sign of (T^ should give reasonable values except in very rare cases. 

The picture for M/M, Fig. is similar. While choosing the second order 
scale is more appropriate than first, causing the second moment to vanish and 
minimizing the third, the zeros and minima of higher moments drift progres- 
sively lower. Such behavior is anticipated by the relatively large value of Aq/q, 
which indicates a wide range of contributing momenta. In this case, higher mo- 
ments are increasingly sensitive to lower q, and the corresponding coefficients 
will progressively increase. We might improve the convergence of the series 
by methods discussed in Sect. ||. For example, choosing q* separately at each 
odd order in as, causing the following even moment to vanish and miminizing 
the subsequent odd moment, with each q* indicating the characteristic scale 
for that moment. An alternative is to resum the entire set of polarization dia- 
grams 0. Regardless, the ability to detect sensitivity to a large range of 
momenta, in addition to the scale itself, by computing the first few logarithmic 
moments is sufficient to warn of large higher order corrections. In this case, it 
suggests using M rather than M in expressions for b and t decays. 

9.2 Quark mass and energy renormalization in lattice NRQCD 

Lattice Nonrelativistic QCD (NRQCD) is an effective field theory designed to 
reproduce the results of continuum QCD for a heavy quark at energies small 
relative to its mass |7|, |7^, |7^ . Higher-dimensional operators provide systematic 
corrections ordered by quark velocity v and lattice spacing a, and account for 
radiative processes above the cutoff, typically around the mass. For a cutoff 
much larger than Aqcd, lattice perturbation theory should give reliable values 
for the coefficients of these operators as well as the renormalization factors 
which connect bare to physical quantities. Ref. demonstrates that this 
expectation is valid, provided one uses a renormalized rather than bare coupling 
constant, and divides link gauge fields by their mean value to remove large 
tadpole contributions peculiar to the lattice. 

Refs. to present calculations of two of these quantities to first order 
in Us', the renormalization factor Zm, which connects the bare lattice heavy 
quark mass to its pole mass, and Eq, the shift from zero of the nonrelativistic 
energy of a heavy quark at rest. Ref. js^, using an action improved to O(w^) 
and O(a^), and to 0{v'^) for spin-dependent interactions, found that first-order 
scale setting produced anomalous results for certain values of the bare mass, 
particularly after tadpole improvement. 

In Tables |] to ^ and Figures ^ to we present new values for the scale for 
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Table 1: Applications of second-order scale setting to several processes. The 
coefficients c„ are defined in Sect. ||. gives the scale set by Eq. (||). 52 gives 
the preferred scale by Eq. ( p^ ) where appropriate, also indicated by boxes. Aq 
measures the range of momentum running through the giuon. 



ci/co 



C2/C0 



12 



Aq 



-logW^ii 
2.448 



3.402/a 



-.691772 0.7076^5 

M/M: 



6.316 0.3194 



-.186421 -.66497 



1.064Vi 



-4.6862 0.09603M 



TjB ^ XueV)/Mf,: 
-4.3163 0.11554M6 

T{t^bW)/M^: 



-5/3 e-5/6Q = 0.4346O 25/9 



17.623 -4.3374 0.27205M 0.38M 



-2.2751 0.32060M^ 5.6848 0.50872 

-5.3382 0.06932Mb 34.410 5.9139 



.0992 -10.531 



23.046 5.3611 



6.0996 -26.477 



0.58534M6 0.35Mb 



-4.2054 0.12213Mf 

r(t ^ bW)/M^t: 
-5.7076 0.05763Mt 

J^'dx[gr{x,Q^)-gr{x,Q')]: 

-2 e-iQ = 0.3679Q 115/18 43/18 

j'dx[Frix,Q')-Frix,Q')]: 



0.75502Mt 



-8/3 e-4/3g = 0.2636g 155/18 3/2 





0.96/a 



0.1 IM^ 



0.084Mh 



0.14Mf 
0.34Mt 
0.28g 
0.16Q 
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Figure 10: The moments |(( [log((7^) — log((7^)]"))|^^" as functions of 
log(g^M2) for n = 1 to 8 (left to right) for T{B X„eFe) over the MS 
mass Mb- The vertical hne indicates the choice for the scale log(g^/M^) us- 
ing Eq. (p^. By choosing the second-order prescription such that the second 
moment vanishes, \og{q^/M^) is either near the minimum or the zero for all 
higher moments, minimizing higher-order terms in Eq. (^5|). 
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Figure 12: The BLM scale q* for the pole mass renormalization factor Zm as a 
function of the bare lattice mass aMp in NRQCD without tadpole improvement. 
The first order solution determines q* for all values. 

a variety of bare quar k masses Mq, both with and without tadpole improvement. 
By applying Eq. ( |l3|) in regions where appropriate, we obtain a reasonable scale 
for all values of Mq, correcting the anomalies observed in Ref. As ex- 

pected, there is a significant reduction in the scale after tadpole improvement. 
The tadpole contributions to these renormalizations are quadratically divergent 
in the inverse lattice spacing, and so are generally large and sensitive to large 
momenta. Tadpole improvement is designed to remove the bulk of these contri- 
butions, and so reduces the typical scale from 2 - 4 to 0.5 - 1.5 in units of the 
inverse lattice spacing a. 

10 Conclusions 

In this paper we have derived a method which incorporates information from 
higher orders into the general prescription of Ref. [Q for choosing the optimal 
scale q* for the strong coupling constant a^. We find that it corrects erroneous 
scales where the leading term or terms are anomalously small. 

The extended prescription states that Eq. (|l3|) determines the optimal scale 
q* when the argument of the square root is positive. When it is not, the first 
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aM, 







Figure 13: The BLM scale q* for the pole mass renormaHzation factor Zm as 
a function of the bare lattice mass oMq in NRQCD with tadpole improvement. 
The first order solution determines q* between oMq = 2.00 and 3.50, the second 
order elsewhere. Circles indicate the appropriate scale; triangles indicate the 
first-order solution in regions where it does not apply. 
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Figure 14: The BLM scale q* for the energy shift Eq as a function of bare 
lattice mass oMo in NRQCD without tadpole improvement. The first order 
solution determines q* for all values. 
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Figure 15: The BLM scale q* for the energy shift Eq as a function of bare 
lattice of the bare lattice mass aMp in NRQCD with tadpole improvement. The 
first order solution determines q* between oMo = 0.80 and 1.70, the second 
order elsewhere. Circles indicate the appropriate scale; triangles indicate the 
first-order solution in regions where it does not apply. 



30 



Table 2: The BLM scale for the pole mass renormalization factor Z,„ for several 
values of the bare lattice mass qMq in NRQCD without tadpole improvement. 
aqi gives the scale set by Eq. p) in units of the inverse lattice spacing. 
gives the preferred scale by Eq. ( |l3|) where appropriate. The parameter n is set 
to ensure the stability of heavy quark propagator evolution in simulations |]76| . 



n 


aMa 


(/) = {Zra - l)/a. 


(/log((ag)2)) 


(flog^iiaqr)) 


aql aq2 


1 


20.00 


0.4679(39) 


1.343(11) 


4.091(62) 


4.202(71) 




17.50 


0.4860(34) 


1.364(10) 


4.196(65) 


4.068(58) 




15.00 


0.5125(30) 


1.3907(85) 


4.211(41) 


3.883(44) 




12.50 


0.5410(46) 


1.426(14) 


4.242(63) 


3.736(65) 




10.00 


0.5880(35) 


1.463(12) 


4.384(68) 


3.469(43) 




7.00 


0.7057(17) 


1.5423(58) 


4.880(32) 


2.982(15) 




5.00 


0.8624(14) 


1.7153(48) 


5.538(29) 


2.7034(87) - 




4.00 


1.0071(27) 


1.9021(70) 


6.062(32) 


2.571(11) 


2 


4.00 


1.0177(23) 


1.9264(80) 


6.164(33) 


2.577(11) 




3.50 


1.1268(23) 


2.0710(83) 


6.614(53) 


2.507(10) 




3.00 


1.2853(21) 


2.2859(68) 


7.351(35) 


2.4333(74) - 




2.70 


1.4119(19) 


2.4540(79) 


7.911(34) 


2.3846(73) - 




2.50 


1.5188(23) 


2.6141(74) 


8.390(52) 


2.3646(66) - 




2.00 


1.9018(22) 


3.1745(79) 


10.214(32) 


2.3039(53) - 




1.70 


2.2751(24) 


3.7546(75) 


12.039(56) 


2.2823(43) 




1.60 


2.4384(24) 


4.0086(79) 


12.857(39) 


2.2750(41) - 




1.50 


2.6320(22) 


4.3141(78) 


13.745(37) 


2.2694(37) - 


3 


1.40 


2.8804(23) 


4.7587(76) 


15.073(37) 


2.2842(34) 




1.20 


3.5010(22) 


5.7986(83) 


18.143(39) 


2.2891(30) - 




1.00 


4.4915(19) 


7.4780(68) 


23.109(31) 


2.2990(19) - 


5 


0.80 


6.3033(34) 


10.720(11) 


32.131(63) 


2.3405(24) 
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Table 3: The BLM scale for the pole mass renormalization factor for several 
values of the bare lattice mass uMq in NRQCD with tadpole improvement, aql 
gives the scale set by Eq. p) i n units of the inverse lattice spacing. gives 
the preferred scale by Eq. ( |l3| ) where appropriate. The parameter n is set to 



ensure the stability of heavy quark propagator evolution in simulations ||76| . 


n 


aMo 


(/) = (Z„, - l)/a. 


(/log((ag)2)) 


(/log^((ag)2)) 


aql 


aq* 


1 


20.00 


-0.2381(39) 


-0.385(11) 


-0.367(62) 


2.246(61) 


1.34(11) 




17.50 


-0.2224(34) 


-0.371(10) 


-0.278(65) 


2.301(60) 


1.240(91) 




15.00 


-0.1996(30) 


-0.3530(85) 


-0.286(41) 


2.422(61) 


1.262(73) 




12.50 


-0.1773(46) 


-0.333(14) 


-0.293(63) 


2.56(12) 


1.29(13) 




10.00 


-0.1416(35) 


-0.323(12) 


-0.223(68) 


3.14(16) 


1.21(12) 




7.00 


-0.0566(17) 


-0.3242(58) 


0.066(32) 


17.6(1.7) 


0.95(13) 




5.00 


0.0386(14) 


-0.3018(48) 


0.335(29) 


0.0201(31) 


0.75(17) 




4.00 


0.1126(27) 


-0.2881(70) 


0.413(32) 


0.278(12) 


0.650(59) 


2 


4.00 


0.1232(23) 


-0.2638(80) 


0.515(33) 


0.343(13) 


0.471(81) 




3.50 


0.1722(23) 


-0.2664(83) 


0.586(53) 


0.461(12) 






3.00 


0.2381(21) 


-0.2783(68) 


0.737(35) 


0.5574(85) 






2.70 


0.2828(19) 


-0.3107(79) 


0.781(34) 


0.5773(84) 






2.50 


0.3180(23) 


-0.3261(74) 


0.807(52) 


0.5988(74) 






2.00 


0.4183(22) 


-0.4581(79) 


0.846(32) 


0.5783(57) 






1.70 


0.4899(24) 


-0.6166(75) 


0.765(56) 


0.5329(44) 


0.57(12) 




1.60 


0.5131(24) 


-0.7058(79) 


0.699(39) 


0.5027(42) 


0.723(23) 




1.50 


0.5376(22) 


-0.8143(78) 


0.518(37) 


0.4689(37) 


0.835(17) 


3 


1.40 


0.5795(23) 


-0.8755(76) 


0.542(37) 


0.4698(34) 


0.839(15) 




1.20 


0.6212(22) 


-1.2529(83) 


-0.043(39) 


0.3648(28) 


1.009(13) 




1.00 


0.6518(19) 


-1.9239(68) 


-1.139(31) 


0.2286(16) 


1.152(11) 


5 


0.80 


0.6964(34) 


-3.009(11) 


-3.277(63) 


0.1153(15) 


1.293(24) 
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Table 4: The BLM scale for the energy shift Eq for several values of the bare 
lattice mass aAfp in NRQCD without tadpole improvement, ag* gives the scale 
set by Eq. (^) in units of the inverse lattice spacing. gives the preferred scale 
by Eq. (|l^) where appropriate. The parameter n is set to ensure the stability 
of heavy quark propagator evolution in simulations |f76| . 



n 


aMa 


(/) = Eo/as 


(/log((ag)2)) 


(/log^((ag)2)) 


aql aq2 


1 


20.00 


2.25711(48) 


2.5490(13) 


10.8100(38) 


1.75884(55) 




17.50 


2.27657(44) 


2.6011(12) 


10.9496(40) 


1.77051(52) - 




15.00 


2.30333(45) 


2.6694(13) 


11.1328(39) 


1.78507(54) - 




12.50 


2.33892(69) 


2.7653(22) 


11.3882(87) 


1.80606(92) - 




10.00 


2.39196(65) 


2.9062(23) 


11.7729(84) 


1.83583(92) - 




7.00 


2.50419(41) 


3.1988(14) 


12.5623(49) 


1.89398(57) - 




5.00 


2.64980(48) 


3.5784(15) 


13.5643(49) 


1.96446(59) - 




4.00 


2.77360(93) 


3.8995(38) 


14.4035(99) 


2.0198(15) 


2 


4.00 


2.77159(78) 


3.8981(28) 


14.380(11) 


2.0203(11) 




3.50 


2.8585(10) 


4.1118(34) 


14.952(11) 


2.0529(13) 




3.00 


2.97012(88) 


4.3973(34) 


15.670(12) 


2.0965(13) 




2.70 


3.0574(11) 


4.6137(38) 


16.228(11) 


2.1266(15) 




2.50 


3.12485(90) 


4.7822(41) 


16.660(12) 


2.1494(15) 




2.00 


3.3492(12) 


5.3483(43) 


18.116(15) 


2.2221(16) 




1.70 


3.5394(13) 


5.8316(43) 


19.302(14) 


2.2791(15) 




1.60 


3.6163(13) 


6.0299(45) 


19.843(16) 


2.3018(16) 




1.50 


3.7057(14) 


6.2486(45) 


20.406(17) 


2.3236(16) 


3 


1.40 


3.7865(16) 


6.4417(61) 


20.873(18) 


2.3411(21) 




1.20 


4.0175(17) 


7.0245(61) 


22.344(17) 


2.3970(20) 




1.00 


4.32658(96) 


7.8282(35) 


24.418(11) 


2.4711(11) 


5 


0.80 


4.6581(20) 


8.7418(74) 


26.784(27) 


2.5557(23) 
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Table 5: The BLM scale for the energy shift Eq for several values of the bare 
lattice mass qMq in NRQCD with tadpole improvement, aql gives the scale set 
by Eq. (|) in units of the inverse lattice spacing, agrj gives the preferred scale 
by Eq. (1^) where appropriate. The parameter n is set to ensure the stability 
of heavy quark propagator evolution in simulations |f76| . 



n 


aMa 


(/) = Eo/as 


(/log((ag)2)) 


(/log^((ag)2)) 


aql 


aq^ 


1 


20.00 


1.05283(48) 


-0.3998(13) 


3.2048(38) 


0.82706(52) 






17.50 


1.04985(44) 


-0.4027(12) 


3.2027(40) 


0.82549(50) 






15.00 


1.04669(45) 


-0.4076(13) 


3.1970(39) 


0.82306(51) 






12.50 


1.04040(69) 


-0.4143(22) 


3.1878(87) 


0.81948(89) 






10.00 


1.03060(65) 


-0.4272(23) 


3.1758(84) 


0.81281(90) 






7.00 


1.00820(41) 


-0.4643(14) 


3.1150(49) 


0.79431(56) 






5.00 


0.97428(48) 


-0.5243(15) 


2.9833(49) 


0.76409(58) 






4.00 


0.94100(93) 


-0.5878(38) 


2.8305(99) 


0.7318(15) 




2 


4.00 


0.93900(78) 


-0.5891(28) 


2.807(11) 


0.7307(11) 






3.50 


0.9137(10) 


-0.6502(34) 


2.671(11) 


0.7006(13) 






3.00 


0.87572(88) 


-0.7311(34) 


2.444(12) 


0.6588(13) 






2.70 


0.8466(11) 


-0.7996(38) 


2.266(11) 


0.6236(15) 






2.50 


0.82102(90) 


-0.8590(41) 


2.111(12) 


0.5927(15) 






2.00 


0.7312(12) 


-1.0621(43) 


1.583(15) 


0.4837(16) 






1.70 


0.6442(13) 


-1.2576(43) 


1.018(14) 


0.3768(15) 


0.7950(59) 




1.60 


0.6056(13) 


-1.3421(45) 


0.830(16) 


0.3302(14) 


0.8460(65) 




1.50 


0.5641(14) 


-1.4439(45) 


0.566(17) 


0.2781(14) 


0.9029(73) 


3 


1.40 


0.4953(16) 


-1.6171(61) 


0.088(18) 


0.1955(16) 


0.986(12) 




1.20 


0.3523(17) 


-1.9501(61) 


-0.802(17) 


0.06282(99) 


1.106(24) 




1.00 


0.13779(96) 


-2.4285(35) 


-2.035(11) 


0.0001489(93) 


1.23(11) 


5 


0.80 


-0.3161(20) 


-3.4380(74) 


-4.629(27) 


230.0(8.4) 


1.416(75) 
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order formula in Eq. (^) applies. The choice of sign for the second-order solution 
should be apparent either from continuity, or by checking that the solution 



minimizes the next higher (cubic) moment in Eq. (55) if it is available. In 
addition, higher moments give a measure of the range Aq of momenta which 
flow through the gluon, and can confirm that the q* chosen in either case is 
indeed typical. Large values for the relative range Aq/q can indicate large 
higher-order contributions even when q* is large. 

Our second-order prescription has several advantages. It requires a simple 
extension to the calculation, either numeric and analytic, needed to implement 
the first-order prescription, requiring only computation of an additional log- 
arithmic moment. Calculation of higher moments can then help to further 
characterize the diagram and confirm the scale choice. It can also identify cases 
where the first two terms are anomalously small, though such cases are ap- 
parently rare. It is appropriate regardless of the number of loops included in 
the running coupling, and is not limited to the strong interactions. Finally, it 
remedies erroneous scales in a variety of processes. 
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